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Abstract 



Thermodynamic stability of statistical systems requires that susceptibili- 
ties be semipositive and finite. Susceptibilities are known to be related to the 
fluctuations of extensive observable quantities. This relation becomes non- 
trivial, when the operator of an observable quantity is represented as a sum 
of operators corresponding to the extensive system parts. The association of 
the dispersions of the partial operator terms with the total dispersion is ana- 
lyzed. A special attention is paid to the dependence of dispersions on the total 
number of particles N in the thermodynamic limit. An operator dispersion 
is called thermodynamically normal, if it is proportional to N at large values 
of the latter. While, if the dispersion is proportional to a higher power of N, 
it is termed thermodynamically anomalous. The following theorem is proved: 
The global dispersion of a composite operator, which is a sum of linearly in- 
dependent self-adjoint terms, is thermodynamically anomalous if and only if 
at least one of the partial dispersions is anomalous, the power of iV in the 
global dispersion being defined by the largest partial dispersion. Conversely, 
the global dispersion is thermodynamically normal if and only if all partial 
dispersions are normal. The application of the theorem is illustrated by sev- 
eral examples of statistical systems. The notion of representative ensembles is 
formulated. The relation between the stability and equivalence of statistical 
ensembles is discussed. 
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I. INTRODUCTION 



Stability of statistical systems and the fluctuations of observable quantities are known to 
be intimately related. The fluctuations can be characterized by the corresponding suscepti- 
bilities, such as specific heat, isothermal compressibility, or longitudinal magnetic suscepti- 
bility. The susceptibilities are connected with the dispersions of the operators representing 
observable quantities. In what follows, we shall deal with the so-called extensive observables, 
whose averages are proportional to the total number of particles N, when N is large [1,2]. 
The existence of the thermodynamic limit is assumed, when N is asymptotically large, such 
that iV — > oo. 

Note that susceptibilities can also be connected with the fluctuations of intensive ther- 
modynamic variables, such as pressure and temperature [3,4]. However, in this paper we 
shall consider only the fluctuations of extensive observables. 

For stable statistical systems in equilibrium, the susceptibilities are positive and finite, 
which follows from their relations to the dispersions of the corresponding operator observ- 
ables [5] or, on the general thermodynamic level, stems from the second law of thermody- 
namics [6]. The susceptibilities may become divergent only at the points of second-order 
phase transitions, which, however, by definition, are the points of instability. Really, at the 
point of a phase transition, one phase becomes unstable, as a consequence, it transforms to 
another, stable, phase. After the phase transition has occurred, all susceptibilities in the 
stable phase go finite. 

The fluctuations of extensive observables, related to the corresponding operator disper- 
sions, can be classified onto two types, according to their dependence on the total number of 
particles N in the given statistical system, when the number N is large, such that N >> 1. 
This implies that the thermodynamic limit is assumed. The fluctuations are called thermo- 
dynamically normal, when the related operator dispersion is proportional to N . Conversely, 
if the operator dispersion is proportional to N a , with a > 1, then the related fluctuations 
are termed thermodynamically anomalous. 

The finiteness of susceptibilities in stable equilibrium systems means that the correspond- 
ing fluctuations are thermodynamically normal. Oppositely, the divergence of susceptibil- 
ities at the critical points shows that the fluctuations of the related extensive observables 
are thermodynamically anomalous. In a stable system, outside phase transition points, all 
susceptibilities are finite, which tells that the fluctuations of all extensive observables are 
thermodynamically normal. 

It is worth warning the reader that thermodynamically normal or anomalous fluctua- 
tions have nothing to do with the normal, that is, Gaussian distributions. Thermodynamic 
normality or anomaly are the notions describing the thermodynamic behaviour of the re- 
lated operator dispersions with respect to the total number of particles. In calculating the 
corresponding averages any quantum or classical probability measures, of arbitrary nature, 
can be employed. 

In the present paper, general relations between the fluctuations of observables and the 
stability of statistical systems are studied. The emphasis is on the case, which is not a 
standard one, when the observable quantities are represented as sums of several terms, cor- 
responding to macroscopic parts of the system. Then the relation between the fluctuations 
of the partial terms and the fluctuations of the global observables is not evident. A general 
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theorem is rigorously proved, connecting the behaviour of fluctuations of global and partial 
observables. This theorem is briefly formulated in the Abstract and its mathematically rigor- 
ous formulation is given in Section III. The direct interrelation between the thermodynamic 
behaviour of fluctuations and stability is emphasized. It is also shown that the stability 
of statistical systems is intricately connected with the notions of symmetry breaking and 
ensemble equivalence. 



II. FLUCTUATIONS OF OBSERVABLES AND STABILITY 

In quantum statistical mechanics, observable quantities are represented by self-adjoint 
operators from the algebra of observables. As is explained in the Introduction, only exten- 
sive observables are considered in the paper. Fluctuations of the observable quantities are 
characterized by the related operator dispersions. Let A be an operator representing an 
extensive observable quantity. Its dispersion is 

A 2 (i) = < A 2 > - < A > 2 , (1) 

where the angle brackets, as usual, denote statistical averaging. 

The dispersions of the operators, representing extensive observables, are directly con- 
nected with the associated susceptibilities, which can be measured. Thus, the fluctuations 
of the Hamiltonian H, quantified by its dispersion A 2 (H), define the specific heat 

where E =< H > is internal energy, N is the total number of particles in the system of 
volume V, and T is temperature. Here and in what follows, the Boltzman constant is set to 
unity, ks = 1. The fluctuations of the number of particles are described by the dispersion 
A 2 (N) of the number-of-particle operator N, yielding the isothermal compressibility 

1 (dV\ A 2 (N) 

in which P is pressure, N ee< N >, and p = N/V is the average particle density. In magnetic 
systems, with the Zeeman interaction — //q Y2% B • Sj of the operator spins Si with an external 
magnetic field B, the fluctuations of the magnetization M a =< M a > are described by the 
dispersion A 2 (M a ) of the magnetization operator M a = fj, J2iLi Sf, which results in the 
longitudinal magnetic susceptibility 

NT K J 

In the notation, used above, /xq = ^7s, with 75 being the gyromagnetic ratio for a particle 
of spin S. In what follows, we shall use the system of units setting to unity the Planck 
constant % = 1. 

The specific heat (2), isothermal compressibility (3), or magnetic susceptibility (4) are 
the examples of the susceptibilities associated with the fluctuations of observables. These 
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thermodynamic characteristics are readily measured in experiments. At the points of phase 
transitions, the susceptibilities can diverge, since such points are the points of instability. 
But for stable equilibrium system, the susceptibilities are always positive and finite for all 
N, including the thermodynamic limit, when N — > oo, V — > oo, so that p = N/V — > const. 
In principle, it is admissible to imagine the situation, when a phase transition occurs not 
merely at a point but in a finite region of a thermodynamic variable [7], inside which region 
the system remains unstable and displays a divergent susceptibility. Such a case, however, 
is quite marginal, and rarely, if ever, happens for real statistical systems. In any event, as 
soon as the phase transition is over, so that the system becomes stable, all susceptibilities 
go finite. 

The following picture summarizes the above consideration. The extensive observables of a 
statistical system are represented by Hermitian operators. The fluctuations of an observable, 
represented by an operator A, are quantified by the operator dispersion A 2 (A), whose ratio 
A 2 (A)/N to the total number of particles characterizes the associated susceptibility. For a 
stable system, the latter must be semipositive and finite, while if it is divergent or negative, 
the system is unstable. This can be formulated as a necessary stability condition 

A 2 (A) 

< —±-L < oo . (5) 

The ratio A 2 (A)/N plays the role of a generalized susceptibility, related to the operator A. 
Examples of condition (5) are the stability conditions on the specific heat (2), isothermal 
compressibility (3), and magnetic susceptibility (4), according to which 

< Cy < OO , < K T < OO , < Xa < OO . (6) 

These thermodynamic characteristics are usually strictly positive at finite temperature, be- 
coming zero only at zero temperature. 

In this way, the dispersion of the operator A, representing an extensive observable, has 
to be proportional to the number of particles: 

A 2 (i) ocN . (7) 

Then the dispersion is called thermodynamically normal. The thermodynamic limit is as- 
sumed here, so that N >> 1. When Eq. (7) is not satisfied, so that A 2 (A) oc N a with a > 1, 
the dispersion is called thermodynamically anomalous. Respectively, the fluctuations of the 
related observable, characterized by the dispersion A 2 (A), are termed thermodynamically 
normal, provided Eq. (7) is valid, and they are named thermodynamically anomalous if Eq. 
(7) does not hold. 

In stable systems, the fluctuations of observables are always normal, and the correspond- 
ing susceptibilities are finite. These susceptibilities can be measured in experiment, either 
directly or through other measurable quantities. For example, the isothermal compressibility 
can be measured through the sound velocity 

2 1 fdP\ 1 
s =-177 = 777777-' ( 8 ) 
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where m is the particle mass. The compressibility can also be found from the central value 
of the structural factor 

S(0) = — 2 = P Tk t . (9) 
ms z 

And the structural factor 

S(k) = l + pJ[g(r)-l]e~^dr, (10) 
in which g(r) is the pair correlation function, can be measured in scattering experiments. 



III. THEOREM ON TOTAL FLUCTUATIONS 

In some cases, the operators of observables have the form of the sum 

a = ea ( n ) 

i 

of self- adjoint terms A\. As has been stressed above, we consider here only extensive ob- 
servables, such that the statistical average < A > is proportional to the total number of 
particles N, when the thermodynamic limit N — > 00 is implied. All parts A, are assumed to 
have the same dimension as A and also to be the operators of extensive observables, so that 
< Ai >oc N. For example, A\—K and A 2 = W could be kinetic and potential energies for 
a system of N particles. Then Eq. (11) would give the Hamiltonian H = K + W. Or one 
can consider the operator of the number of particles N = N + Ni as a sum (11) composed 
of the operators of condensed particles, Nq, and of noncondensed particles, Aq, for a system 
with Bose-Einstein condensate. For each of the terms, one may consider partial fluctuations 
quantified by the dispersions A 2 (Aj). Then of the principal interest is the problem how 
the partial dispersions A 2 (A;) are correlated with the total dispersion A 2 (A)? For instance, 
could it be that some of the partial dispersions are thermodynamically anomalous, while 
the total dispersion remains thermodynamically normal, so that the system as a total stays 
stable? The answer to such questions is given by the following theorem. 

Theorem. Let the operator A of an extensive observable quantity be represented as 
a sum of linearly independent self-adjoint operators Ai, being of the same dimension and 
also representing extensive observables, such that < Aj >oc N in the thermodynamic limit. 
Then the global dispersion A 2 (A) is thermodynamically anomalous, so that A 2 (A) oc N a 
with a > 1, if and only if at least one of the partial dispersions A 2 (Aj) is thermodynamically 
anomalous. The power a in the dependence A 2 (A) oc N a , as N — > 00, is defined by the 
largest power of all partial dispersions A 2 (Aj). Conversely, the global dispersion A 2 (A) is 
thermodynamically normal, such that A 2 (A) oc A in the thermodynamic limit, if and only 
if all partial dispersions A 2 (Aj) are thermodynamically normal. 

Proof. First, let us note that it is meaningful to consider only linearly independent 
terms in the sum (11), since in the opposite case, when some of the terms are linearly 
dependent, it is straightforward to express one of them through the others, so that to reduce 
the number of terms in sum (11). For concreteness, in the following proof, the representatives 
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of observables are called operators, which assumes the case of a quantum system. Of course, 
the same argumentation is valid for classical systems as well, for which one just has to 
replace the term "operator" by the term "classical random variable". 
The dispersion for the operator sum (11) can be written as 

A 2 (i) = £ A 2 (i.) +2^cov(iU,) , (12) 

i i<j 

where the covariance 

cov(A, Aj) = ^< AAj + AjAi > - < Ai>< Aj > (13) 

is employed. The latter enjoys the symmetry property 

cov(Ai,Aj) = cov(Aj,Ai) . 

The dispersions are, by definition, semipositive, while the covariances can be positive as well 
as negative. 

It is sufficient to prove the theorem for the sum of two operators, when 

A 2 (i, + A,) = A 2 (A) + A 2 (AA + 2cov(i,, A,) . (14) 

This follows from the simple fact that any sum of terms more than two can always be 
redefined as a sum of two new terms. We assume that in Eq. (14), where i ^ j, both terms 
are operators but not classical functions. If one of the terms were just a classical function, 
then we would have a trivial equality 

A 2 (A; + const) = A 2 (Ai) , 

with the left-hand and right-hand sides being simultaneously either thermodynamically nor- 
mal or anomalous. 
The elements 

Gij = cov(A, Aj) , (15) 

having the properties an = A 2 (Ai) > and = <jji, form the covariance matrix [cr^]. This 
matrix is symmetric. For a set of arbitrary real- valued numbers Xi, with i = 1,2, ... ,n, 
where n is an integer, one has 



n - 

< Yl {Ai- < Ai >)xi 



>= (TijXiXj > . (16) 



The right-hand side of equality (16) is a semipositive quadratic form. The theory of quadratic 
forms [8] tells us that a quadratic form is semipositive if and only if all principal minors of its 
coefficient matrix are non-negative. Thus, the sequential principal minors of the covariance 
matrix [cry], with i,j = 1,2, ... ,n, are all non-negative. In particular, 



a ii a jj ~ a ij a ji — 
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This, because of the symmetry = a^, takes the form 



Hence, the correlation coefficient 



possesses the property 



A« = (17) 



A?-<1- 



The equality A?- = 1 holds true if and only if A\ and Aj are linearly dependent. The 
sufficient condition is evident, since if Aj = a + bA i: with a and b being any real numbers, 
then cry = 6<Tjj and cr^- = ftVjj, thence A^ = b/\b\, from where A?- = 1. To prove the 
necessary condition, let us assume that A^- = 1. Therefore Ay = ±1. Let us consider the 
dispersion 

f 1 ±4^) = 2(l±A y )>0. 



The value Ay = 1 is possible then and only then, when 
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A 



The dispersion can be zero if and only if 

= const , 

which implies that the operators Ai and Aj are linearly dependent. In the same way, the 
value Xij = — 1 is possible if and only if 

A 2 (-^= + -^1 = . 



And this is admissible then and only then, when 



A.? A. j 

H — = const , 



which again means the linear dependence of the operators Ai and Aj. As far as these 
operators, by assumption, are linearly independent, one has 



K < 1 • (18) 



This inequality is equivalent to 



a J < auajj , 
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which, employing notation (15), becomes 

cov(A t ,A 3 )\ 2 < A 2 (A)A 2 (i,-)- (19) 

Now, equality (14) can be represented as 

A 2 (Ai + Aj) = an + (jjj + ■l\, lX '(T,,(T !j , (20) 

where, as is shown above, < 1. Altogether there can occur no more than four following 
cases. First, both partial dispersions an = A 2 {Ai) and ajj = A 2 {Aj) are normal, so that 
an oc iV and a^ oc N. Then, from Eq. (20) it is obvious that the total dispersion A 2 {Ai + 
Aj) oc iV is also normal. Second, one of the partial dispersions, say an oc N, is normal, 
but another one is anomalous, a^ oc N a , with a > 1. From Eq. (20), using the inequality 
(1 + a)/2 < a, one has A 2 (Ai + Aj) oc N a . That is, the total dispersion is anomalous, 
with the same power a as ajj. Third, both partial dispersions are anomalous, such that 
an oc N ai and ajj oc iV aj with different powers, say 1 < ctj < aj. Then Eq. (20), with 
taking account of the inequality (c^ + <X/)/2 < aj, shows that A 2 (Ai + Aj) oc . Hence, 
the total dispersion is also anomalous, with the power aj of the largest partial dispersion ajj. 
Fourth, both partial dispersions are anomalous, an oc c 2 N a and ajj oc c 2 N a , where q > 
and Cj > 0, with the same power a. In that case, Eq. (20) yields A 2 (Ai + Aj) = CijN a , 
where 

Cij oc - Cj) 2 + 2ciCj(l + \j) > , 

which is strictly positive in view of inequality (18). That is, the total dispersion is anomalous, 
having the same power a of iV as both partial dispersions. After listing all admissible 
cases, we see that the total dispersion is anomalous if and only if at least one of its partial 
dispersions is anomalous, with the power of iV of the total dispersion being equal to the 
largest power of partial dispersions. Conversely, the total dispersion is normal if and only if 
all its partial dispersions are normal. This concludes the proof of the theorem. 

This theorem was, first, announced, without proof, in Ref. [9]. The proof, presented 
above, is rather general, being valid for arbitrary operators and statistical systems. The 
theorem can be applied to any system. For instance, this can be a multicomponent system, 
where the index % in Eq. (11) enumerates the components. In recent years, much attention is 
given to systems with Bose-Einstein condensate (see review articles [10-12]). The problem of 
fluctuations in such systems has received a great deal of attention, with a number of papers 
claiming the existence of anomalous fluctuations everywhere below the condensation point 
(see discussion in Ref. [13]). In the following sections, the examples of Bose-condensed 
systems will be considered. In addition to being naturally separated into the condensed 
and noncondensed parts, Bose systems can also display the coexistence of several coherent 
topological modes [14-23]. Another possibility is the coexistence of atoms in several internal 
states, which, e.g., has been studied in collective Raman scattering [24]. 



IV. IDEAL BOSE GAS 

The uniform ideal Bose gas below the condensation temperature is known to exhibit 
anomalous number-of-particle fluctuations [25,26]. Here, this case will be briefly recalled for 
the purpose of illustrating the above theorem. 
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The condensation temperature of the ideal uniform Bose gas is 

2/3 



nn 27T 

T r = — 



m 



P 



(21) 



_C(3/2). 

where C(3/2) = 2.612. Below this temperature, the number-of-particle operator is the sum 

N = N + iVx (22) 
of the terms corresponding to condensed and noncondensed particles, respectively, 

N = aja , N 1 = Y^ a \ a k , 

where a k and are the creation and annihilation operators of Bose particles with momentum 
k. 

The dispersion for the total number-of-particle operator N can be calculated by means 
of the derivative over the chemical potential //, so that 



A 2 (N) = T 



dN 



(/i - -0) . 



The average number of particles N =< N > is given by the sum 

N = N + N 1 

of condensed, 

N = <ala >= (e'^-iy 1 , 

and noncondensed, 



N EE < N > = 



N 



93/2 



(23) 



(24) 



(25) 



(26) 



particles, where \i — > — 0, 



/ 2tt 
mT 



and the Bose-Einstein function is 



. , 1 r°° zu n ~ l , 

9n{z) = / du . 

1 in) Jo 



e"- — z 



Let us stress that the terms N and Ni in the sum (23) are linearly independent. Differen- 
tiating the sum (24), one has the total dispersion 



A 2 (N) = A 2 (N ) + A 2 (N 1 ) , 



(27) 
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with the partial dispersions 



From Eqs. (25) and (26), we find the dispersion for condensed particles, 

A 2 (iV ) = N (l + N ) , (28) 

and for noncondensed particles, 

A 2 (iV x ) = -1L g 1/2 (e?») , (29) 

where fi — > —0. As far as the existence of Bose-Einstein condensate presupposes that the 
number of condensed particles N is macroscopic, that is, proportional to N, then from 
Eq. (28) and the relation N oc TV > 1, we have A 2 (iV ) oc N 2 . Expression (29) in the 
thermodynamic limit possesses an infrared divergence caused by the integral 

. . l r°° du 

#1/2(1) oc 



■ - *J Um.irt. ^ ^ 

in which 

k 2 1 
2mT L 

with L oc V 1 / 3 . Consequently, #1/2(1) oc L/X T . Thus, dispersion (29) diverges at finite 
temperatures as 

A 2 (A>0 oc (mT) V 4/3 . (30) 

In this way, both dispersions for the number-of-particle operators of condensed as well as 
noncondensed particles are anomalous: 

A 2 (iV ) oc N 2 , A 2 ^) oc iV 4/3 . 

As a result, the total dispersion (27) is also anomalous, A 2 (AQ oc iV 2 , with the power of iV 
given by A 2 (iVo). 

The anomalous dispersion A 2 (A) leads, according to Eq. (3), to the divergence of the 
isothermal compressibility, as kt oc N, everywhere below T c , except T = 0. But the system 
with a divergent compressibility is not stable. Therefore, the ideal uniform Bose gas below 
the condensation temperature (21) is a pathological object, being unstable in the whole 
region < T < T c . In other words, such a gas does not exist as a stable statistical system 
[13]- 

It is worth emphasizing that the anomalous fluctuations of the condensate can be cured 
by breaking gauge symmetry as will be explained below. However the fluctuations of noncon- 
densed particles remain anomalous, with the dispersion A 2 (A^!) oc N^ 3 in both ensembles, 
grand canonical as well as canonical [25,26]. Therefore, the instability of the ideal uniform 
Bose gas below T c is not an artifact caused by the choice of an ensemble, but a property 
peculiar to this system. 
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V. INTERACTING BOSE GAS 



There exists a popular myth that the number-of-particle fluctuations of noncondensed 
particles in an interacting Bose gas below T c remain anomalous, corresponding to the dis- 
persion A 2 (A r 1 ) oc iV 4 / 3 , of the same type as that for the ideal Bose gas (see discussion in 
Ref. [13]). If this were true, then according to the theorem of Section 3, the total dis- 
persion A 2 (N) would also be anomalous, with the power of iV not smaller than 4/3. This 
would imply that the isothermal compressibility diverges at least as k t oc iV 1 / 3 . Hence the 
system as a whole would be unstable. In turn, this would mean that there are no stable 
statistical systems with Bose-Einstein condensate. Such a conclusion, of course, would be 
too radical, because of which it is necessary to reconsider the procedure of calculating the 
number-of-particle dispersions for Bose-condensed systems. 

Let us consider a weakly interacting Bose gas at low temperatures, when the Bogolubov 
theory [27-29] is applicable. The main points of this theory are as follows. One starts with 
the standard Hamiltonian 

H = J V> f (r) (-^--Jj ^(r) dr + i J ^ t (r)^ t (r')$(r - r')^(r')V(r) drdr' (31) 

in terms of the Bose field operators ip(r) and '(r). The interaction potential is assumed to 
be symmetric, such that $(— r) = $(r), and soft, allowing for the Fourier transformation 

$ ( r ) = 77 £ ^ e * k ' r > = / $ (r)e- lk ' r dr . 

v k J 

The condensate is separated by means of the Bogolubov shift 

^(r)=^ + ^(r), (32) 

in which 

= ~% , ^i(r) = £ °fcV?fc(r) , (33) 

and, keeping in mind a uniform system, the expansion is over the plane waves <p k (r) — 
e tk ' r /^/V. The gauge symmetry of Hamiltonian (31) is broken by setting a = \^Nq. As- 
suming that N Q pa N, one omits from the total Hamiltonian the terms of the third and fourth 
order with respect to the operators a k of noncondensed particles, where k ^ 0. Retaining 
only the terms up to the second order in a k , one gets the quadratic Hamiltonian 

H 2 = ]- Np<5> + ]T uj k a\a k - fiN + ]- ]T A k (a\a)_ k + a_ k a k ) , (34) 

1 fc^O 1 k^O 

in which the notation for the quantities 

k 2 

uJ k = — + p($ + $fc) — A* (35) 
2m 

and 
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A k = p$ fe (36) 

is employed. 

The quadratic Hamiltonian (34) is diagonalized by means of the Bogolubov canonical 
transformation 

a k = u k b k + v*_ k tf_ k , 

in which 

2 2 1 

fc 2efc 2efe ' " 2e k 2e k 

and £fc is the Bogolubov spectrum 

e k = y/<4 - A\ . (37) 

The condensate separation through the Bogolubov shift (32) is meaningful only when the 
particle spectrum (37) touches zero at k — 0, which gives 

fi = p$ • (38) 

Thus, one comes to the Bogolubov Hamiltonian 

H B = E + Y,£ k blb k -iiN , (39) 

fc/0 



with the ground-state energy 



E = \ Np% - \ YJ^k ~ e k ) ■ (40) 

Using the chemical potential (38), for the spectrum (35) one has 

k 2 

u k = — + p$ k . (41) 
2m 

With the diagonal Bogolubov Hamiltonian (39), it is easy to find the normal, 

n k = < a\a k > , (42) 

and anomalous, 

a k =< a k a^ k > , (43) 

averages. We have 

n k = ^-(l + 27r k )- ± (44) 

and 
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a k = -^(l + 27r k ), (45) 

where 

7lk = < b\b k > = (e^-\y l . (46) 

Now let us turn to investigating the number-of-particle fluctuations. In the Bogolubov 
approximation, the number-of-particle operators for condensed, N , and noncondensed, N±, 
particles are uncorrelated, so that 



Hence, their covariance 
Therefore 



< iVoiVi > = < N >< Nx > . (47) 
cov(iVo,iVi) = . 

A 2 (N) = A 2 (N ) + A 2 (iYx) . (48) 



Calculating the dispersion A 2 (iVi) for the number-of-particle operator of noncondensed 
particles 

Nl = J2 a l a k > 

one has to work out the four-operator expression < a\a k a q a q > or, after involving the Bo- 
golubov canonical transformation, one needs to treat the four-operator terms < b k b k b q b q >. 
Such four-operator products are reorganized by means of the Wick decoupling, which yields 

A 2 M = E { (l + 2 -^f) 7r fc (l + n k ) + ^} . (49) 

Here the notation 

[p®k~ 

Ck = \ 

V m 

for the effective sound velocity is used, which enters the Bogolubov spectrum (37) as 



k 2 x 2 



t kk)2 + \^) ' (50) 



Replacing in Eq. (49) the summation by integration, one gets an infrared divergence of 
the type iV / dk/k 2 . Limiting here the integration by minimal k min = 1/L, with L oc N 1 ^ 3 , 
one gets A 2 (iV"i) oc N 4 ^ 3 , which is anomalous. Remaining in the frame of the discrete wave 
vectors k does not save the situation, and the dispersion A 2 (iVi) stands anomalous. But, as 
follows from the theorem of Sec. Ill, the anomalous partial dispersion yields the anomalous 
total dispersion A 2 (iV), which in the present case is evident from Eq. (48). As a result, the 
compressibility (3) diverges as k t oc iV 1//3 , which implies the instability of the system as a 
whole. Thus one would come to the strange conclusion that stable Bose-condensed systems 
do not exist. 
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However, the conclusion on the appearance of anomalous fluctuations in Bose systems, 
derived from Eq. (49), is not correct. The mistake here is in the following. A basic point of 
the Bogolubov theory is the contraction of the total Hamiltonian (31) to the quadratic form 
(34), omitting all terms of the order higher than two with respect to the operators a k of 
noncondensed particles. The Bogolubov theory is a second- order theory with respect to a k - 
Being in the frame of a second-order theory imposes the restriction of keeping only the terms 
of up to the second order when calculating any physical quantities, and omitting all higher 
order terms. In working out the dispersion A 2 (A r 1 ), one meets the fourth-order terms with 
respect to a k . Such fourth-order terms are not defined in the second-order approximation. 
The calculation of the fourth-order expressions in the second-order approximation is not 
self-consistent, i.e., it is incorrect. 

A correct calculation of A 2 (N) in the frame of the Bogolubov theory can be accomplished 
in the following way. By invoking the relations (3), (9), and (10), we have 

A 2 (iV) = iv|l+p|[( ? (r)-l]rfr} . (51) 

The pair correlation function is 

g(r 12 ) = \< ^ t (ri)V t (r 2 )V(r 2 )V'(ri) > , (52) 
P 

where r 12 = i"i — r 2 . 

For the field operators, one assumes the Bogolubov shift (32), which taking into account 
that in the thermodynamic limit the condensate operator ip becomes a classical number, 
can be written as 

V>(r) = r? + ^i (r) , (53) 
where the first term is the Bogolubov order parameter 

v = < ^( r ) > = < > , (54) 



which can be set as rj = y/po, with p = N /V. Here rj does not depend on r for a uniform 
system under consideration. 

The pair correlation function (52) can be simplified by invoking the Wick decoupling. 
This, however, must be handled with care. A delicate point is that the Wick decoupling and 
the Bogolubov shift (53) do not commute with each other. In the present context, the Wick 
decoupling is equivalent to the Hartree-Fock-Bogolubov approximation. The latter does not 
commute with the Bogolubov shift. Thus, accomplishing, first, the Bogolubov shift in the 
pair correlation function (52), and then using the Hartree-Fock-Bogolubov approximation 
for the operators ^i(r), or, what is the same, the Wick decoupling for the operators a k , with 
k 7^ 0, we obtain 

g(r 12 ) = 1 + Re[/>i(ri,r 2 ) +<7i(ri,r 2 )] + \ [|pi(n, r 2 )| 2 + K(n, r 2 )| 2 j . (55) 
p A p z 1 J 

Here the Hartree-Fock-Bogolubov approximation for ipi(r) is employed, resulting in 

< Vl(ri)Vi(ri)Vi(r 2 ) > = 0, 
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because of the condition < ipi(r) >= 0, and in 

< V t (ri)V'i(r2)V'i(r2)V'i(ri) > = p\ + |pi(ri, r 2 )| 2 + K(n, r 2 )| 2 . 
The notation is used for the normal average 

pi(ri, r 2 ) = < ^(r 2 )^i(ri) > (56) 
and for the anomalous average 

(7i (ri, r 2 ) =< ^i(r 2 )^i(ri) > (57) 
in the real space. These averages are related, by means of the Fourier transforms 

Mr., r 2 ) = / ^ , Mr, r 2 ) = / ^ , 

with the normal and anomalous averages (42) and (43), respectively, in the momentum 
space. 

Note that function (55) possesses the correct limiting behaviour 

lim gr(r 12 ) = 1 . 

But, if one, first, would make the Hartree-Fock-Bogolubov approximation for the operators 
ip(r) and, after this, would substitute the Bogolubov shift (53), then one would get another 
correlation function with a wrong limiting behaviour, as is explained in the Appendix A. This 
is because the usage of the Wick decoupling, and Hartree-Fock-Bogolubov approximation, 
for the operators, represented as sums of several terms, is correct if and only if all terms 
in the sum possess the same commutation relations. However, in the Bogolubov shift (53), 
the field operators ip(r) and ipi(r) do have the same Bose commutation relations, but the 
term 77 does not enjoy such relations. Consequently, the proper way of action is to realize, 
first, the Bogolubov shift (53) and only after this to invoke the Hartree-Fock-Bogolubov 
approximation for the operators ipi(r). The inverse order, as is explained in the Appendix 
A, is not correct. 

For the pair correlation function (55), we find 

2 Po v 1 \ , 1 f ( 2 , 2\ 



J [g(r) - 1] dr = hm(n k + a k ) + - J (n\ + a 2 k ) 



In the frame of the Bogolubov theory, we have to set po = P and to omit the terms of the 
order higher than two with respect to the operators a k of noncondensed particles. This means 
that the terms n\ and a\ are to be omitted. Therefore, the number-of-particle dispersion 
(51) in the Bogolubov theory is 



A 2 (iV) = N 
Employing Eqs. (44) to (46), we get 



1 + 2 lim(n k + a k ) 



(58) 



limK + a k ) = \{^-l 
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where 

k^o " v m 

with 



c = lim Cfe 



$ = lim = / $(r) dr . 

— ^0 J 

Then dispersion (58) becomes 



which is, of course, normal, as it should be for a stable system. Respectively, the isothermal 
compressibility 



A 2 (N) 

pTN pmc 2 



k t = = — (60) 



is finite. 

According to the theorem of Sec. Ill, if the total dispersion (59) is normal, then both 
dispersions of the number-of-particle operators for condensed, A 2 (A r ), as well as for non- 
condensed, A 2 (Ni), particles must be normal. Anomalous fluctuations can arise solely as a 
result of wrong calculations, when, e.g., one considers the fourth-order terms n\ and a\ in 
the second-order Bogolubov theory. 



VI. SYSTEMS WITH CONTINUOUS SYMMETRY 

It is easy to show that the same fictitious anomalous fluctuations appear, not only for 
Bose systems, but for arbitrary systems, when one treats the Hamiltonian in the second-order 
approximation, but intends to calculate fourth-order expressions. This immediately follows 
from the analysis of susceptibilities for arbitrary systems with continuous symmetry, as has 
been done by Patashinsky and Pokrovsky [30]. Following Ref. [230], one may consider an 
operator A = A(<p), which is a functional of a field ip. Let this operator be represented as a 
sum A = A +A 1 , where the first term is quadratic in the field if, so that A oc </?V> while the 
second term depends on the field fluctuations Sip as A 1 oc Sip^Sip. Let the system Hamiltonian 
be taken in the hydrodynamic approximation, where only the terms quadratic in the field 
fluctuations 5ip are retained. The dispersion A 2 (A) oc Nx is proportional to a longitudinal 
susceptibility x- The latter is given by the integral / C(r)dr over the correlation function 
C(r) = g(r) — 1, with g(r) being the pair correlation function. Calculating A 2 (A), one meets 
the fourth-order term < 5ip^5ip8ip^8ip >. For the quadratic hydrodynamic Hamiltonian, such 
fourth-order terms are decoupled by resorting to the Wick theorem. Then one finds 

<?(') « ^i) ( 61 ) 
for any dimensionality d > 2. Consequently, 

X oc J C(r) dr oc iV^ 2 )/ 3 
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for 2 < d < 4. Hence the dispersion is 

A 2 (A) oc N X oc jV (d+1)/3 . (62) 

For d — 3, this gives A 2 (A) oc iV 4//3 , that is, the same anomalous dispersion as A 2 (iV) for 
Bose systems. But this implies that the related susceptibility diverges as x oc iV 1 / 3 , which 
tells that the considered system is unstable. If this would be correct, it would mean that 
there are no stable systems with continuous symmetry. For instance, there could not exist 
magnetic systems, described by the Heisenberg Hamiltonian. Liquid helium also could not 
exist stable system. 

The existence or absence of anomalous fluctuations does not depend on the statistical 
ensemble used. Thus, in the frame of the same calculational procedure, the particle fluctua- 
tions are the same, being either anomalous or normal, depending on the chosen procedure, 
for all ensembles, whether canonical, grand canonical, or microcanonical [31]. 

It is worth emphasizing that such fictitious anomalous fluctuations arise not just at a 
phase transition point, which would not be surprising, but everywhere below this point, 
in the whole region of existence of the considered system. That is, everywhere below the 
phase transition points such systems would not be stable. As is evident, such a strange 
conclusion is physically unreasonable. Fortunately, the explanation for the occurrence of 
anomalous fluctuations is rather simple: They arise solely due to an incorrect calculational 
procedure, when the fourth-order terms are treated by a second-order theory, such as the 
hydrodynamic approximation. No anomalous fluctuations happen, if all calculations are 
done self-consistently, being defined in the frame of the given approximation. 

Another popular way of incorrectly obtaining thermodynamically anomalous particle 
fluctuations for systems with continuous symmetry is as follows. One uses the representation 

^(r) = e^ (r) yn^) (63) 

for the field operator, in which n(r) = ^(r)^(r) is the operator of particle density and </?(r) 
is the phase operator. The latter is assumed to be Hermitian in order to preserve the correct 
definition of the density operator, 

^ f (r)^(r) = yn(r) e -^ +(r)+ ^ (r) yn(r) = n(r) . 

It is easy to show that from the representation (63) it follows that the density and phase 
operators are canonically conjugated, satisfying the commutation relation 

[n(r), t/3(r')] = iS(r — r') . 

For the first-order correlation function, one has 

< ^(r)^(0) > = < ^Jn{r)n(0) exp {-i [<p(r) - 0(0)]} > . 

Then one assumes that the temperature is asymptotically low, T — > 0, such that there are no 
density fluctuations, and one can replace the operator n(r) by its average p(r) =< n(r) >. 
This is equivalent to the usage, instead of the representation (63), of the representation 

V>(r) = y^)e^. (64) 
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One also supposes that the phase fluctuations are very small, so that one can employ the 
following averaging: 

< exp {-i [y3(r) - £(0)]} > = exp j- ± < [y3(r) - y3(0)f >} . (65) 

As a result, the first-order correlation function reduces to 

< V f (r)V(0) > = P(r) exp {- i < [y3(r) - y?(0)] 2 >} . 

Treating ip(r) as a small quantity, one also expands the exponentials in powers of <p(r). 
Similarly, one treats the second-order correlation functions. Finally, one comes to the same 
expressions as in Eqs. (61) and (62), with the thermodynamically anomalous fluctuations, 
A 2 (Ai) oc A 4 / 3 , for the three-dimensional space. 

The main mistake in such calculations is the same as has been made above. All calcula- 
tions have been based on the assumption that both the density a and phase fluctuations are 
rather weak, so that the hydrodynamic approximation could be invoked. The latter implies 
that all statistical averages are treated in the hydrodynamic approximation, with a Hamil- 
tonian quadratic in the operators. For instance, it is well known [32] that Eq. (65) is valid 
solely for quadratic Hamiltonians. For finding A 2 (Aq), one needs to consider the fourth- 
order terms in phase operators. Of course, there is no sense in calculating the forth-order 
terms in the frame of a second-order theory, such as the hydrodynamic approximation. 

Moreover, the representations (63) and (64), as such, are principally incorrect. This 
is shown in the Appendix B. A correct definition of the phase operator requires a much 
more elaborate technique, as can be inferred from the review articles [33-36]. Since the 
representations (63) and (64), actually, do not exist, all conclusions derived on their basis, 
even involving no further approximations, are not reliable. 



VII. BREAKING OF GAUGE SYMMETRY 

In Section IV, considering the ideal uniform Bose gas, we found that its particle fluctu- 
ations are thermodynamically anomalous, with the corresponding dispersion A 2 (A) oc A 2 . 
This anomaly is due to the condensate fluctuations, since A 2 (A ) oc A 2 . Really, for an ideal 
uniform gas, one has 

A 2 (A)=5> fc (l+n fc ). (66) 

k 

From here, separating the terms with k = and k ^ 0, we get 

A 2 (A ) = A (l + N ) , A 2 (A X ) = n k (l + n k ) . 

Since A oc A, we find A 2 (A ) oc A 2 . 

The situation can be made even more dramatic by generalizing it to the case of inter- 
acting particles. To this end, let us consider an interacting system that can be treated 
by perturbation theory starting with a mean-field approximation, such as the Hartree-Fock 
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approximation. In the frame of the latter, the particle dispersion can be shown [37] to have 
the same form as in Eq. (66). Then, irrespectively of the concrete expression for the mo- 
mentum distribution of particles nk, the global dispersion A 2 (TV) will be thermodynamically 
anomalous because of the anomalous term A 2 (iV ) oc iV 2 . Hence, one could conclude that 
all systems with the Bose-Einstein condensate would be unstable. 

One often states that the appearance of this anomaly is the defect of the grand canonical 
ensemble. However, this is not correct. As is mentioned in Section IV, the anomalous 
condensate fluctuations are fictitious and can be removed by breaking the gauge symmetry. 

Hohenberg and Martin [38] noticed that the appearance of such fictitious divergences 
is a common feature of theories possessing gauge symmetry, but breaking the latter would 
eliminate the divergences resulting from the condensate fluctuations. Ter Haar [25] showed 
explicitly how the anomalous condensate fluctuations can be removed after breaking the 
gauge symmetry for an ideal uniform Bose gas. In the present section, we demonstrate that, 
in general, the gauge-symmetry breaking eliminates the anomalous condensate fluctuations 
for arbitrary systems, whether interacting or not. 

A known method for lifting a system symmetry of any nature is the method of infinites- 
imal sources, introduced by Bogolubov [29,39]. There are also several other methods of 
symmetry breaking, as is reviewed in Ref. [40]. In the case of gauge symmetry, one has to 
be cautious by chosing the way of its breaking. The standard method of infinitesimal sources 
may not always lead to the desired symmetry breaking, as is shown by a counterexample in 
the Appendix C. 

To break the gauge symmetry in a Bose system, one has to resort to the Bogolubov shift 
[29,39]. The latter, keeping in mind the most general statistical system, whether equilibrium 
or nonequilibrium, uniform or nonuniform, writes as 

i>(r,t) = r)(r,t)+if> 1 (r,t) , (67) 

where t is time. The first term here is the condensate wave function, assumed to be not iden- 
tically zero in the presense of the Bose-Einstein condensate. The second term in Eq. (67) is 
the field operator of noncondensed particles, satisfying the same Bose commutation relations 
as ip(r,t). The correct separation of condensed and noncondensed particles presupposes the 
orthogonality condition 

Jr ] *(r,t)if; 1 (r,t)dr = , (68) 

which exculdes the double counting of the degrees of freedom. In what follows, just for 
brevity, we shall write ^(r) instead of ip(r, t), understanding that, generally, the time variable 
t does enter the dependence of the field operator, ^(r) = ip(r,t). 

For the theory of Bose systems, it is extremely important to specify the spaces of states, 
which the field operators are defined on. Thus, the field operators if>(r) and if>'(r) are defined 
on the Fock space F{if;) generated by the operator if; '(t). This means the following [41]. 
There exists a vacuum state |0 >, for which 

ip(r)\0 > = . (69) 

The Fock space T{if;) is the space of all states 

oo -y . n 

v=Yl ~n / /»( r i> • • • ' r «) II ^( r i) dr i\° > > 
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in which / n (ri, • • • , r n ) is a square- integrable function symmetric with respect to the permu- 
tation of any pair of its variables. 

It is easy to notice that the state |0 >, which is a vacuum state for ^(r), is not a vacuum 
for ip±(r), since 

V>i(r)|0>= -77(r)|0 > ^ 0. 
Consequently, there should exist another state |0 >i satisfying the condition 



^i(r)|0>i=0, 



(70) 



being a vacuum for ip\(r). In turn, the state |0 >i, which is a vacuum for ipi(r), is not a 
vacuum for ^(r), as far as 

^(r)|0>i=77(r)|0>i^O. 

The Bogolubov shift (67) is a particular case of canonical transformations [42]. The 
operators ^(r) and ipi(r) can be connected with each other by means of the transformation 



C = exp < / rf* (r)tp(r) — rj(r)ip'(r) 



dr 



and its inverse 



C- 1 = exp 



r)*(r)ip(r) — r}(r)ip'(r) 



dr 



Using these transformations, one has 
and 

^i(r) = C-^{v)C . 
Then it becomes clear that the vacuum for ipiir) is 

|0 > l = c-^o > . 



(71) 



(72) 



(73) 



(74) 



(75) 



The vacua |0 > and |0 >i are mutually orthogonal. This can be shown by employing 
the Baker-Hausdorff formula, which for two operators A and B, whose commutator [A, B\ 
is proportional to the unity operator, reads as 



e A+B = e A e B exp 



A, B 



Using this for transformation (72), we have 

C- 1 = exp || 77(r)^ f (r) drj exp j- J ?f (r)^(r) drj expj-- J \r](r)\ 2 rfrj . (76) 
Acting on the vacuum |0 >, we find 

C _1 |0 > = ex p{-^ J\v(r)\ 2 drj exp [J rj{r)^ (t) dr ] 1° > • ( 77 ) 
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This is nothing but the coherent state [43] , being the eigenstate of the destruction operator, 

V>(r)|?7>= 77(r)|?7> , (78) 
and having in the coordinate representation [44] the form 

\v > — Vo ex P^J i](r)^(r) drj |0 > , (79) 



with the normalization factor 

\Vo\ 

Respectively, the condensate wave function 

r)(r) = < r]\ip(r)\r] > 

is nothing but the coherent field related to the coherent state \i] >. 
In this way, the vacuum (75) is the coherent state (79), 

|0 >! = CT^O > = 1 77 > . (80) 

The scalar product of the vacua |0 > and |0 >i is 

< 0|0 >i = < 0|?7 > = exp j-^ J \r](r)\ 2 dr j . (81) 

By its definition, the condensate wave function gives the condensate density 

Po(r) = Mr)| 2 . (82) 

The number of condensed particles 

N = [ p (r) dr , (83) 



in the presence of the condensate, is not zero, but is macroscopic in the sense that iVo oc 
N — > oo. Therefore the scalar product 

< 0|0 >! = exp (- I N ) (84) 



becomes zero in the thermodynamic limit, 

< 0|0 >i ~ (AT — > oo) . (85) 

This tells that the vacua |0 > and |0 >i are asymptotically orthogonal. The Fock spaces 
and generated from the related vacua, are orthogonal to each other, except just 

the sole state |0 >i= 1 77 >, which is the vacuum for J-^i) and the coherent state, defined 
by Eq. (78), in F^ip). However, having the sole common state for two infinite-dimensional 
spaces means the intersection of zero measure. Moreover, the influence of this intersection 
is eliminated by means of the orthogonality condition (68). 
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Thus, there are two different vacua |0 > and |0 >i and two mutually orthogonal Fock 
spaces and J-{ipi)^ generated by the field operators tp^ and respectively. The 

operator (7f) transforms ^(ipi) into F^ip), while the operator (72) transforms Fijp) into 
!F{ip-i). There is no self-adjoint operator C + that would be defined on the same space as 
C. Therefore the operator C is nonunitary and the transformations (73) and (74) cannot 
be treated as unitary. The field operators ip and ip\ are defined on different spaces. One 
says that such operators realize unitary nonequivalent operator representations of canonical 
commutation relations [45]. 

Breaking the gauge symmetry by the Bogolubov shift (67), one, actually, passes from 
the Fock space J 7 ^) to the space !F{ipi). Since the left-hand and right-hand sides of Eq. 
(67) are defined on different spaces, this equation should be understood as a transformation 

ip(r) — ► rj(r) + ^i(r) . 

Separating the zero-momentum mode for a uniform Bose gas, with replacing this term by a 
nonoperator quantity, 

/ a ° i— 
Vo = -j= -»■ VPo , 

as has been done in Section V, is mathematically equivalent to the Bogolubov shift [46]. 
The representation of the operators of observables, expressed through the field operators ipi, 
and defined on the Fock space can be called the Bogolubov representation. 

In the Bogolubov representation, the operator of condensed particles, according to Eqs. 
(82) and (83), is a nonoperator quantity, N = N . Hence, the dispersion of the latter is 
zero, A 2 (iV ) = 0. Consequently, the dispersion of the total number-of-particle operator 

A 2 (N) = A 2 (7V0 

is completely defined by the dispersion of the operator Ni of noncondensed particles. Thus, 
the anomalous iV 2 dispersion of the condensate particles is removed in the Bogolubov rep- 
resentation. 

Considering the ideal uniform Bose gas of Section IV in the Bogolubov representation, we 
do not meet the iV 2 -anomalous condensate fluctuations. Nevertheless, particle fluctuations, 
characterized by the dispersion A 2 (iVi) oc iV 4//3 , remain thermodynamically anomalous. 
That is, this gas, anyway, is unstable. This conclusion does not depend on whether the 
grand canonical or canonical ensemble has been used. Of course, in the latter, where the 
total number of particles is fixed, the related dispersion is not defined. However, one can 
calculate the compressibility 

'd 2 F\~ l 




, x dV 2 



i 



where F is free energy. For the ideal uniform Bose gas below T c , one has [2] dP/dV = 0, 
hence, kt — > oo, which implies instability. The latter is an intrinsic feature of the uniform 
ideal Bose gas [13]. Including particle interactions stabilizes the gas, as is shown in Section V. 
The ideal Bose gas can also be stabilized by trapping it in an external confining potential, 
such as the harmonic potential [47,48], though not all power-law potentials are able to 
stabilize the system [49]. 



22 



The message of this section is that accurately defining the symmetry properties of the 
given system helps to avoid the appearance of unphysical instabilities. Although there also 
exist systems, such as the ideal uniform Bose-condensed gas, which are intrinsically unstable. 

VIII. NOTION OF REPRESENTATIVE ENSEMBLES 

The consideration of the previous Section VII demonstrates the importance of accurately 
defining the system under investigation. It is not sufficient to chose a statistical ensemble, 
but often it is also necessary to formulate additional conditions specifying the features of 
the given system, thus, avoiding the appearance of spurious instabilities. For instance, one 
can take the grand canonical ensemble without breaking the gauge symmetry or one may 
employ the grand canonical ensemble with the gauge symmetry breaking. This means that, 
in general, there may exist not just the sole grand canonical ensemble or the sole canonical 
one, but there can exist several such ensembles. This problem of the ensemble nonuniqueness 
is just another way of formulating the problem of the nonuniqueness of the Fock space and 
of the existence of unitary nonequivalent operator representations, which is explained in the 
previous Section VII. 

Thus, for the correct description of a physical system, it is necessary to equip the chosen 
statistical ensemble by additional conditions required for accurately taking account of the 
system features. Only such an equipped ensemble will correctly represent the considered 
system, that is, will be a representative ensemble. 

The idea of the representative ensembles goes back to Gibbs himself [50], who mentioned 
the necessity of taking into account all additional information known about the considered 
system, such as the system symmetry, the existence of integrals of motion, and so on. The 
importance of employing representative ensembles for an adequate description of statistical 
systems was emphasized by ter Haar [51,52]. A detailed discussion of mathematical tech- 
niques, required for the correct definition of representative ensembles, can be found in the 
review papers [40,53]. In the language of reduced density matrices, the latter have to satisfy 
specific constraints in order to correctly represent a given statistical system [54]. 

Systems, exhibiting Bose-Einstein condensation, serve as a very good example demon- 
strating the importance of taking into account their specific features in order to correctly 
describe their behaviour. Rich properties of these systems require to be very attentive in for- 
mulating the corresponding representative ensemble. Forgetting to impose the appropriate 
constraints, specifying the system properties, may lead to self-inconsistent calculations and 
the appearance of spurious instabilities. In Section V, the example was given of a weakly- 
interacting equilibrium uniform Bose gas. Now we shall formulate a general approach to 
Bose systems with arbitrarily strong interactions, being, in general, nonuniform and not 
necessarily equilibrium. We shall stress the constraints that are compulsory for defining a 
self-consistent theory, which, for equilibrium systems, results in a representative ensemble, 
free of fictitious instabilities. 

First of all, as is explained in Section VII, we have to break the gauge symmetry by 
means of the Bogolubov shift, replacing the field operator ip(r,t), acting in the Fock space 
F(ip), by the operator 

4>(r,t)=r)(r,t) + Mr,t) , (86) 
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defined on the Fock space In what follows, we shall again omit the time variable 

in order to simplify the notation. The first term in the right-hand side of Eq. (86) is the 
condensate wave function and the second term is the field operator of noncondensed parti- 
cles. The replacement tp(r) — > ^(r) yields to the passage from the operator representation 
on the Fock space to the unitary nonequivalent operator representation, the Bogol- 

ubov representation, on the space only if the condensate wave function i](r,t) is not 

identically zero. 

The energy operator has now to be expressed through the field operators (86), which 
yields the Hamiltonian 

H = J ft(r) ^ + Uj ^(r) dr+ 

+ i J ^ t (r)^ t (r / )$(r - r')^(r')^(r) drdr' , (87) 
in which U — U(r, t) is an external field. The corresponding Lagrangian is 

L = J — ip(r) dr - H . (88) 

It is important to stress that, contrary to a system without condensate, where there is 
just one field operator variable vp, in a Bose-condensed system, there appear two variables 
i] and ipi, or one can take as two variables i] and ip. The condensate wave function defines 
the condensate density (82). The operator of the total number of particles 

N = J ft(r)$(r) dr (89) 

is expressed through ip. Respectively, there are two normalization conditions. One condition 
is for the condensate wave function normalized to the number of condensed particles 



iVr 



= jW)?dv. (90) 



And another normalization condition is for ip normalized to the total number of particles 
N =< N >, i.e., 

N = J < ^(r)^(r) > dr . (91) 

Here and everywhere in this section, the angle brackets imply the averaging over the Fock 
space Fiipi). 

Hamiltonian (87), with the field operator (86), contains the terms linear in ipi, because 
of which the average < ipi > may be nonzero. However, a nonzero < ip\ > would, in general, 
lead to the nonconservation of quantum numbers, such as spin and momentum, which would 
be unphysical. Therefore, it is necessary to impose the constraint for the conservation of 
quantum numbers, 

< ^(r) > = . (92) 
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In this way, three conditions are to be valid for a Bose-condensed system, two normalization 
conditions (90) and (91), and the quantum- number conservation constraint (92). 

The most general procedure of deriving the equations of motion is by looking at the 
extrema of the action, under the given additional conditions. In our case, the effective 
action is 



A[r], Vi] = J (1 + fi N + fiN + A) dt . 



(93) 



Here, L is the Lagrangian (88). The second and third terms in the integral (93) preserve 
the normalization conditions (90) and (91). And the role of the term 



A = y [A(r)v4(r) + AWi(r) 



dr 



(94) 



is to satisfy the quantum- number conservation constraint (92). The Lagrange multipliers 
A(r) have to be chosen so that to cancel in Eq. (87) the terms linear in The absence 
of such linear terms in the Hamiltonian, as is known [42] , is necessary and sufficient for the 
validity of condition (92). By introducing the effective grand Hamiltonian 



H[ V , Vi] = H-tioN -tJiN-A 
and the resulting Lagrangian 

L[V, H = J 
the effective action (93) can be rewritten as 

A[ v , Vi] = J L[ v , Vi] dt 



d d 
V*(r)i ^ rj(r) + ip\{r)i — ^i(r) 



dr - H[r}, ^i] , 



(95) 



(96) 



(97) 



According to the standard prescription, the equations of motion are obtained from the vari- 
ational principle determining the extremum of the action functional (97). These variational 
equations are 



5A[r), ip 



5r]*(r, t) 

where, for generality, the time variable is written explicitly, and 



(98) 



. 



(99) 



From Eqs. (95), (96), and (97), it follows that Eqs. (98) and (99) are identical to the 
variational equations 



1 dt^ r,t) ~ SrTfat) ' 



(100) 



with the effective grand Hamiltonian (95), and 
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. d 5H[tjM 

(5^1 (M) 



(101) 



Explicitly, Eq. (100) is 



■ 9 , s 



2m 



+ U -e)r](r) + 



+ J $(r-r') [\r,{v')\ 2 r,{T)+X{v,v')\ dv' , 



(102) 



where e = /i + p and again, for short, the time dependence is omitted. Equation (101) 
yields 

d / V 2 \ 

+ J $(r-r') [|7 ? (rO| 2 ^i(r)+r ? *(rO^(r)V'i(rO+^(rO^(r)V'i t (rO+X(r,r^ dr' . (103) 
Here the notation 

X(r,r') = VlWitr'Mr) + ^(r>(r>i(r)+ 
+7 7 *(r / )Vi(r / )V'i(r) + ^(r>i(r>i(r) (104) 
is used. Averaging Eq. (102), we obtain the equation for the condensate wave function 

' V 2 ^ 



*^(M) 



2m 



+ U-e\ r](r) + 



+ J $(r - r') [p(t')ti(t) + Pl (r, r')^(r') + ^(r, r')r l *(r')+ < ^(r>i(r>i(r) 



> 



dv' , 
(105) 

in which the total density of particles 

p(r) = p (r) +pi(r) 

is the sum of the condensate density (82) and of the density of noncondensed particles 

pi(r) = < ^(r)Vi(r) > ; 
also the notation is used for the normal density matrix 

Pl (r,r') = <^(r')^i(r)>, 

and the so-called anomalous density matrix 

<7i(r,r') = < ^!(r')^i(r) > , 
which is nonzero because of the broken gauge symmetry. 
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It is not our goal to study here particular consequencies of the approach sketched above. 
The sole aim of the example of this section is to illustrate the way of constructing a rep- 
resentative ensemble for a rather nontrivial system. This is done by accurately specifying 
the basic system properties, such as the broken gauge symmetry, normalization conditions, 
and the quantum-number conservation condition. Following the most general procedure 
of action variation, under the specified conditions, one automatically obtains an effective 
Hamiltonian and the related exact equations of motion. It is possible to show [37] that the 
latter guarantee the correct behaviour for the spectrum of collective excitations, the validity 
of all conservation laws, and the absence of unphysical instabilities. 

It may happen in some lower-order approximations that there is no need to invoke all of 
the conditions discussed above. This, for instance, occurs in the Bogolubov approximation 
of Section IV. In this approximation, one assumes that Nq — > N, hence //o — > 0. Also, 
for a uniform gas, the Hamiltonian term of the first order in ipi vanishes itself, while the 
terms of the third and fourth order in tpi are neglected in the Bogolubov second-order 
approximation. Because of this, there is no necessity of introducing the term (94). However, 
all these conditions are to be taken into account when going to higher-order approximations. 
In the other case, the defined ensemble may occur to be nonrepresentative, which can result 
in physical inconsistences and fictitious instabilities. 

Correctly defining a representative ensemble is also crucially important for the problem 
of equivalence of statistical ensembles, which is discussed in the next section. 

IX. PROBLEM OF ENSEMBLE EQUIVALENCE 

The examples of the previous sections show that the stability properties of a system 
can be different in different ensembles. More general, the same physical quantity may be 
different, being calculated in two different ensembles. Does this mean the failure of the basic 
principle of statistical mechanics, stating the equivalence of ensembles for large systems? 
This question is analyzed in the present section. 

First of all, let us stress that, as is clear from the previous sections, a physical system 
and a describing it ensemble do not exist separately, but they are intimately connected. A 
correct formulation of an ensemble does presuppose that it includes the information on the 
main system features. An ensemble, which is adequate for the given physical system, is 
only that, which properly represents the system, that is, a representative ensemble. But if 
there are two representative ensembles for the same system, then, by their definition, they 
must yield identical results for the same physical quantities. In the other case, at least 
one of these ensembles does not correctly describe the system, hence, is not representative. 
Also, in the case of equilibrium, it is meaningful to talk only about stable systems, as far 
as an unstable system cannot be in absolute equilibrium. Thus, in terms of representative 
ensembles, the following statement is straightforward: Two ensembles are equivalent if and 
only if both of them are representative for the given stable system. Conversely, when two 
ensembles are not equivalent, then at least one of them is not representative. An ensemble 
that is not representative for the given system may be representative for some other system. 
However, there is no any reason to require that two ensembles applied to two different 
physical systems be equivalent. Ensemble nonequivalence, vaguely formulated, is a rather 
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artificial nonphysical problem caused by an improper usage of ensembles not representing 
the considered system. 

To be more correct, let us recall that, generally, one distinguishes two types of ensem- 
ble equivalence, thermodynamic and statistical. In thermodynamics, a physical system is 
characterized by thermodynamic potentials, each of which is a function of its natural ther- 
modynamic variables [1-7]. The system is stable, when thermodynamic potentials enjoy the 
property of convexity or concavity with respect to the appropriate variables. The thermody- 
namic potentials, expressed through different thermodynamic variables, are connected with 
each other by Legendre transforms [1-7]. All thermodynamic characteristics are defined as 
derivatives of thermodynamic potentials. When the latter are connected by Legendre trans- 
forms and correspond to a stable (in the sense of the convexity or concavity property of the 
potentials) system, then the thermodynamic characteristics, calculated in different ensem- 
bles, coincide with each other. Summarizing, the concept of thermodynamic equivalence can 
be formulated as follows: 

Thermodynamic equivalence. Two ensembles , representing a stable physical system, are 
thermodynamically equivalent if and only if their thermodynamic potentials are mutually 
connected by Legendre transforms. 

A rigorous proof of this statement for the case of the macrocanonical and canonical 
ensembles can be found in Refs. [55,56]. Several examples of systems with long-range 
interactions have been considered, whose microcanonical entropy is not a concave function of 
energy [55-58] . The internal energy of such systems, though being nonadditive, can be made 
extensive by means of the Kac-Uhlenbeck-Hemmer normalization [59] yielding a well defined 
thermodynamic limit. The canonical free energy is a concave function of inverse temperature, 
but the microcanonical entropy is not a concave function of energy. This does not allow to use 
the Legendre transform in both directions [55,56]. The nonconcavity of the microcanonical 
entropy results in the appearance, for some range of energies, of negative specific heat, while 
in the canonical ensemble specific heat is always positive. Because of this, one tells that, for 
such models with long-range interactions, the microcanonical and canonical ensembles are 
not equivalent. However, a microcanonical ensemble with a nonconcave entropy does not 
represent a stable physical system, i.e., this ensemble is not representative. As is explained 
above, there is no sense to compare nonrepresentative ensembles, which are not obliged to be 
equivalent. To make the microcanonical ensemble representative, it must be complimented 
by the concavity construction rendering stability again. After this, it becomes representative 
and completely equivalent to the canonical ensemble. 

Nonconcave microcanonical entropy and negative specific heat are also known for grav- 
itating systems, as is reviewed in Refs. [60,61]. To avoid the negative specific heat, one 
can again invoke a concavity construction or to use the canonical ensemble. However, con- 
trary to other models with long-range interactions, the energy of gravitating systems, being 
proportional to iV 5//3 , cannot be made extensive, which does not allow the existence of the 
thermodynamic limit. For gravitating systems, the condition of global equilibrium [62] 

— > const < (106) 

is not valid. Therefore, they may be in principle unstable, which makes questionable the 
application for their description of equilibrium statistical mechanics. 
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The notion of statistical equivalence of ensembles is based on the comparison of the 
averages of observable quantities calculated in different ensembles. To concretize this, let 
us consider the operators of observables A defined on a Fock space T . The set of all these 
operators forms the algebra of observables A = {A}. The statistical state is defined [44,63] as 
the set<„4>={<A>} composed of all statistical averages for the algebra of observables. 
The calculation of the averages is defined in the standard way as the trace of A, with a 
statistical operator corresponding to the chosen ensemble. Let us define as < A > M the 
statistical state related to the grand canonical ensemble, with a chemical potential p. For 
short, the dependence of the state on temperature T and volume V is not shown explicitly. 
For instance, the average density is 

P=y, N=<N>, . (107) 

Suppose, we wish to compare the grand canonical and canonical ensembles. Recall that 
the general structure of the Fock space is a direct sum 

T = ®™ =0 H n (108) 

of the n-particle Hilbert spaces TC n . The pertinent mathematical details can be found in 
Refs. [41,42,44,63]. Define a restriction of the operator A on 7i n as A n . Then the statistical 
state in the canonical ensemble can be denoted as < An > p , with a fixed density p and the 
number of particles N. In view of the structure (108), the states < A >^ and < An > p are 
related through the integral 

< A > Mp) = / K(p,x) < A N(x) > p dx , (109) 

J 

in which p = p(p) is a solution of Eq. (107) and N(x) = xV. The kernel K(p,x) is called 
the Kac density. The corresponding states coincide, when in the thermodynamic limit 

K(p,x) -> 5(p-x) . 

Then one has 

< A > Mp) = < A N > p , (110) 

which signifies the statistical equivalence of grand canonical and canonical ensembles. 

Comparing the statistical states, one has to be very cautious, remembering that it may 
happen that there is not just the sole canonical or grand canonical ensemble, but there could 
be several such ensembles depending on additional constraints specifying the properties of 
the considered system. This is related to the nonuniqueness of the Fock space (108) and 
the existence of nonequivalent operator representations, as is discussed in Sections VII and 
VIII. Therefore, one has, first of all, to define the appropriate representative ensembles and 
only after this one can compare the related averages. If at least one of the ensembles is not 
representative, then there is no sense to compare the averages and equality (110) does not 
need to be valid. 

As an example, let us take a Bose-condensed system, which, according to the previous 
sections, can be considered either using an operator representation on the gauge-symmetric 
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space J-{^) or employing the Bogolubov representation on the space J-(tpi), with broken 
gauge symmetry. In the former case, some fictitious instabilities may arise and Eq. (110) 
may become invalid. However, this would not imply nonequivalence of the ensembles, but 
would simply mean that nonrepresentative ensembles are involved. 

Recall as well that a representative ensemble is assumed to represent a stable system. 
For unstable models, Eq. (110) does not have to be always valid. For instance, if we consider 
the ideal Bose gas in a box, which, as has been explained above, is not stable, then there is 
no reason to require that Eq. (110) be true. This is really so below the condensation point 
[64,65], where the Bose-condensed gas becomes unstable. This instability is manifested by 
thermodynamically anomalous density fluctuations. The ideal Bose gas is also shown [65] to 
be unstable with respect to boundary conditions, whose slight variation leads to a dramatic 
change of the spatial particle distribution, even in the thermodynamic limit. This is contrary 
to the behaviour of realistic stable systems, for which the influence of boundary conditions 
disappears in the thermodynamic limit. Changing, for the ideal Bose gas, the boundary 
conditions from repulsive to attractive [65] transforms the Bose-Einstein condensation from 
the bulk phenomenon to a strange surface effect, when the condensate is localized in a narrow 
domain in the vicinity of the system surface, being mainly concentrated at the corners of 
an infinite box. It is clear that a system, in which the condensate is localized somewhere at 
the corners of an infinite volume, is a rather unphysical object. 

Thus, formally comparing two ensembles, one sometimes can arrive at their seeming 
nonequivalence. This, however, in no way invalidates the basic principle of statistical me- 
chanics stating the ensemble equivalence. This just means that at least one of the compared 
ensembles is not representative, which also includes that the system may be intrinsically un- 
stable. The principle of equivalence holds only for representative ensembles, which represent 
stable systems. 

X. CONCLUSION 

The analysis is given of the relation between the stability properties of statistical systems 
and the fluctuations of observable quantities. The emphasis is made on the composite 
observables that are represented by the sums of several terms. The main result of the 
paper is the theorem connecting the global fluctuations of an observable with the partial 
fluctuations of its components. The theorem is general, being formulated for an arbitrary 
operator represented as a sum of linearly independent self-adjoint operators. These operators 
can be associated with the total and partial observable quantities of a statistical system. The 
theorem tells that: The total dispersion of an operator, being a sum of linearly independent 
self-adjoint operators, is thermodynamically anomalous if and only if at least one of the 
partial dispersions is anomalous, with the power of N in the total dispersion defined by the 
largest partial dispersion. Conversely, the total dispersion is thermodynamically normal if 
and only if all partial dispersions are normal. 

The theorem allows us to understand the relation between the fluctuations of partial ob- 
servables and the fluctuations of the total observable. Respectively, the character of partial 
fluctuations turns out to be directly related to the stability of statistical systems. Several 
examples illustrate the practicality of the theorem, helping to avoid wrong conclusions that 
could happen when studying the behaviour of partial observables. In particular, the fluc- 
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tuations of condensed, as well as noncondensed particles, in a Bose-condensed system must 
be normal, if the system is assumed to be stable. In the same way, fluctuations is systems 
with continuous symmetry are also thermodynamically normal. 

Breaking of gauge symmetry helps to eliminate fictitious instabilities arising in Bose- 
condensed systems. Generally, it is crucially important that a system be characterized by 
its representative ensemble. This makes it possible to avoid artificial contradictions in the 
theory and the related unphysical instabilities. One of the basic principles of statistical 
mechanics, the principle of ensembles equivalence, holds only for representative ensembles 
correctly representing stable statistical systems. 
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Appendix A. Noncommutativity of Bogolubov Shift 

This Appendix illustrates the noncommutativity of the Bogolubov shift and the Hartree- 
Fock-Bogolubov approximation (HFB approximation). When one accomplishes in function 
(52), first, the Bogolubov shift (53) and then the HFB approximation for ipi(r), one gets 
expression (55), with the correct limiting behaviour. But in the other way round, employing, 
first, the HFB approximation for ip(r) and, after this, substituting the Bogolubov shift (53), 
one gets 

g(r 12 ) = 1 + ^2 + ^ RefpiO^ra) +<j 1 (r 1 ,r 2 )] + \ {| Pl (n, r 2 )| 2 + K(n, r 2 )| 2 } . 

p pz pA y ) 

The limiting behaviour of this pair correlation function is not correct, since here 

hm g(r 12 ) = 1 + M , 

ri2->oo p<! 

which would be true only when p = 0. But when p 7^ 0) we confront the problem of the 
condensate overcounting. Thence, these procedures are not commutable. And one has, first, 
to introduce the Bogolubov shift (53) and only after this to resort to the HFB approximation. 
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Appendix B. Nonexistence of Phase Operator 

To show that the representation (63) does not exist, we may use the method of reduction 
to absurdity. Suppose that this representation is correct. Then, from the commutation 
relation 

[n(r), 0(r>)} = iS(r - r') , 
we obtain for the number-of-particle operator 



fn{v) 



N = I n(r) dr 

the commutaton relation 



N, <p(r) 



% . 



From here, taking the matrix element with respect to the number basis {\n >}, for which 
N\n >= n\n >, we find 

(n — n') < n\ip(r)\n' > = i5 nn i . 

Setting here n = n', we get the senseless equality i = 0. Thus, the representation (63) does 
not exist. 

Now, suppose that the representation (64) is correct. Then for the density operator, we 
have 

n(r) = ^(r)tl)(r) = p(r) . 
Hence, the number-of-particle operator becomes identical to the total number of particles, 

N = J p(r) dr = N . 

At the same time, there is an exact relation 

>(r),JV] =V«. 

Using this for N = N, we get the senseless equality ^(r) = 0. Hence, the representation 
(64) is wrong. 

In this way, neither representation (63) nor representation (64) are correct. The phase 
operator, defined through these representations, does not exist. To introduce correctly a 
kind of a quasi-phase operator, one should employ the Pegg-Barnett technique [36]. 
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Appendix C. Gauge-Symmetry Breaking 

The simple method of infinitesimal sources may not always break gauge symmetry. To 
illustrate this, it is sufficient to give at least one counterexample. For this purpose, let us 
consider the Hamiltonian 

H = J ^(r)u(r)^(r) dr , 

with a positive function u(r) > 0. This Hamiltonian is invariant under the gauge transfor- 
mation 

iKr) — > e ! >(r) , 

where a is any real-valued number. Hence < ip(r) >= 0. To break the gauge symmetry, 
following the standard method of infinitesimal sources, one adds to the Hamiltonian H a 
term lifting the symmetry. For instance, the Hamiltonian 



H £ = H-e 



A*(r)V>(r) + A(r)^(r) 



dr 



where A(r) is a complex-valued function, is not gauge invariant. The latter Hamiltonian can 
be diagonalized by means of the canonical transformation 

in which the new field operator ip(r) enjoys the same commutation relations as ip(r). Then 
we have 

H £ = E £ + J$(t)u(t)$(t) dr, 

with the notation 

J uj(r) 

For the diagonal in ijj(r) Hamiltonian H £ , one has < ^(r) >= 0. Therefore 

A(r) 



< ip(r) > = e 



uo r 



According to the method of infinitesimal sources, after calculating the averages, one should 
set e — > 0. But then 

< V(r) > -> (e -> 0) , 

because of which the gauge symmetry has not been broken. Contrary to this, the Bogolubov 
shift (67) is a sufficient condition for gauge-symmetry breaking. 
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